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Abstract. Transiting planets provide a unique opportunity to study the atmospheres of extra- 
solar planets. Radiative hydrodynamical models of the atmosphere provide a crucial link between 
the physical characteristics of the atmosphere and the observed properties. Here I present results 
from 3D simulations which solve the full Navier-Stokes equations coupled to a flux-limited diffu- 
sion treatment of radiation transfer for planets with 1,3, and 7 day periods. Variations in opacity 
amongst models leads to a variation in the temperature differential across the planet, while at- 
i— I ' mospheric dynamics becomes much more variable at longer orbital periods. I also present 3D 
radiative simulations illustrating the importance of distinguishing between optical and infrared 
opacities. 
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JL 1. Introduction 

The abundance of close-in planets coupled with the wealth of information that can 
be deduced from the fraction that transit their host stars is now allowing us to probe 
our understanding of the physics of atmospheres through radiative hydrodynamical mod- 
elling. A number of approaches to this problem exist within the community. The physics 
of radiation and hydrodynamics are both extremely important, and there are several 
approaches currently utilized. Below I separate the approaches to both the dynamical 
and radiation aspects. Although this is not an exhaustive list, I believe it encapsulates 
l^s the most widely utilized approaches in the field. 

Dynamical Methods 

f^*i ■ Equivalent Barotropic Equations (2D) 

00 ! Shallow Water Equations (3D) 

O ■ Full Navier-Stokes Equations (3D) 

J> , Radiative Transfer Methods 

J^j ■ Relaxation Methods 

r\ Flux-limited Radiative Diffusion 

C^ ' Coupled Radiation/Thermal Energy Equations 

Full Wavelength Dependent Radiative Transfer 

Ideally, we would solve the full 3D Navier-Stokes coupled with wavelength dependent 
radiative transfer. However, such a task is well beyond our computational resources, so 
groups working in this area select a method from one or both lists with which to tackle 
the problem. Although selecting methods with many built-in assumptions allows for high 
resolution studies, one runs the risk of missing fundamental physics crucial for a full un- 
derstanding of the atmospheres. On the flip side, as with all numerical studies, the more 
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exact your method, the more difficult and longer it takes to simulate. For a comparison 
of methods please see Dobbs-Dixon & Lin (2008) or Showman et al (2008). 

In the following proceedings I will describe results from simulations utilizing the full 
3D Navier-Stokes equations (described in Section (2)). In Section (3), I couple these 
with the flux-limited radiative transfer in simulations which solve a single equation for 
the energy of the radiation and gas components. These simulations exhibit a range of 
behavior for different opacities and orbital periods. I show that the local value of opacity 
is crucial in determining the efficiency of heat re-distribution to the night-side. In Section 
(4) I present results of simulations again solving the full Navier-Stokes equations, but 
now solving separate energy equations for the radiation and gas components. Decoupling 
these equations allows for the introduction of multiple opacities corresponding to both 
incident stellar and local radiation. As I show, this distinction is crucial in reproducing 
the observed upper atmosphere temperature inversions. I conclude in Section (5) with a 
discussion of the results. 



2. Navier-Stokes Equations 

Planetary atmospheres are simulated using a three-dimensional radiative hydrodynam- 
ical model in spherical coordinates (r, <j), 9) . The equations of continuity and momentum 
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where Q is the rotation frequency and the $ = is the gravitational potential. 

As mentioned above, I present two separate sets of simulations in this proceeding. For 
those presented in Section (3), the radiative and thermal portions are combined into a 
single equation. The energy equation can then be written as 
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= -PVu-VF (2.3) 



where e = c v pT is the internal energy density, T is the local gas temperature, c v is the 
specific heat, and F is the radiative flux. The numerical techniques utilized in solving 
Equations (|2TTj) - ([2~3f are described in Dobbs-Dixon & Lin (2008). 

The final ingredient for solving the energy balance in the atmosphere is a closure 
relation linking the flux F back to the radiation energy density. Here we utilize the 
flux-limited diffusion (FLD) approximation of Levermore & Pomraning (1981), where 

F = -A y -V£ fl . (2.4) 

In the above Er — aT^ and A is a spatial and temporally variable flux-limiter providing 
the closure relationship between flux and radiation energy density. The functional form 
of A is given by 

A= 2 + R 9 , (2.5) 



where 



R - PKR (T,P) E R ■ (2 ' 6) 
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FLD allows for the simultaneous study of optically thin and optically thick gas, correctly 
reproducing the limiting behavior of the radiation at both extremes. In the optically 
thick limit the flux becomes the standard 

4acT3 -VT. (2.7) 
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For the optically thin streaming limit, Equation ()2.4|) becomes |F| = cEr. Between 
these limits Equation (J2.5P is chosen to approximate the full radiative transfer models of 
Levermore & Pomraning (1981). 



3. Orbital Period and Opacity Variations 

In this section I present simulations solving Equations (|2.ip - (|2.3[) in a planet similar to 
HD 209458b, but for a variety of opacities and orbital periods. Opacity is a crucial fac- 
tor in determining the efficiency of energy transport through the atmosphere. Therefore, 
I present separate results utilizing both the higher interstellar opacities (Pollack et al. 
(1985) for lower temperatures coupled with Alexander & Ferguson (1994) for higher tem- 
peratures) and also lower line opacities of Freedman et al. (2008). Both sets of opacities 
are temperature and pressure dependent. Although it is now thought that the opacities 
of Freedman et al. (2008) more closely represent the actual planetary opacities, it is clear 
that there is great diversity among the observed properties. Opacity may likely be a 
major factor in this diversity, so understanding how the dynamics and energy transfer 
change with varying opacity is quite important. 

In Figure (pj I show the effect of opacity on the temperature at the photosphere of a 
synchronously rotating planet at an orbital distance of 0.04 AU. As is clear from Figure 
([TJ, the higher the opacity the lower the night-side temperature. Night-side temperatures 
in a planet with high interstellar opacities are around 300 K, while the night-side of a 
planet with lower opacity is around 700 K. This difference can be easily understood 
by considering that the depth at which radiation is deposited changes with opacity. 
For large opacities the stellar energy is deposited high in the atmosphere where cooling 
timescales are relatively short. The flow is not able to effectively advect energy around 
the planet before it cools. For the lower opacity model, the stellar energy is deposited 
deeper in the atmosphere where cooling timescales become comparable to the crossing 
timescale. The efficiency of energy advection is also evident in the lobes that extend 
onto the night-side at both the equator (in the eastward direction) and the poles (in 
the westward direction) . These regions correspond to the largest eastward and westward 
flow velocities respectively, thus the shortest crossing timescales. In general the night-side 
temperature can be calculated by equating the crossing and radiative timescales (Burkert 
et al. (2005)). This gives a opacity dependent temperature that can be expressed as 
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T n = s — . (3.1) 

\ 3TTKdO-Rp ) 

This is explored in greater detail in Dobbs-Dixon & Lin (2008). 

Given the rapidly developing pace of the field it is also worthwhile to examine the 
atmospheric dynamics of planets for a range orbital periods. Given the short timescale 
for synchronization I assume that the rotation and orbital frequencies are tidally locked in 
each model. The problem of non-synchronous planets (whose subsolar points move across 
the planet with time) is also an interesting and relevant question and will be examined in 
Dobbs-Dixon et al. (2008). In Figure J5J I show the results of simulations with orbital and 
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Figure 1. Temperature variations at the photosphere of irradiated planets synchronously ro- 
tating at orbital periods of 3-days. The left panel shows results using large interstellar opacities, 
while the right panel shows results using lower Freedman et al. (2008) opacities. The subsolar 
point is located at zero degrees. Photospheric velocities are shown with the over-plotted arrows. 




Figure 2. Temperature variations at the photosphere of planets with orbital periods of 1 (left 
panel) and 7 (right panel) days. Clear structures associated with stable jets are evident in the 
1-day simulation, while the night-side temperature of the highly variable 7 day simulation is 
much more uniform. 

—4/3 

rotation periods of 1 and 7 days, with incident fluxes scaling as P orh ■ These simulations 
were completed using the opacities of Freedman et al. (2008) and should be compared to 
the P or b — 3 day simulation shown on the right side of Figure |T|) . 

The clearly evident banded structure in the simulations of 1 and 3 days is a direct 
result of the latitudinal Coriolis force —2H, p v c j 1 sin(d). The formation of the equatorial jet 
comes about as eastward moving material feels a force funnelling it toward the equator, 
while the oppositely moving westward flow is pushed toward the poles. Faster rotation 
results in stronger funnelling and narrower jets. 

The behavior of the simulation with an orbital period of 7 days is somewhat different. In 
general the structure of an eastward equatorial jet and westward jets at higher and lower 
latitudes' persists. However, the day-night temperature difference in this simulations is 
smaller. The jet velocity, driven by the temperature differential, is substantial smaller 
and the resulting temperature and velocity structure on the night-side is quite variable. 



4. Decoupled Energy Equations: Temperature Inversions 

As spectral observations become available for an increasing number of transiting plan- 
ets, the presence of temperature inversions in several systems have become evident 
([eg.]Burrows et al. (2007)). Although suggested by Hubeny et al. (2003), prior to these 
observations this notion had not been pursued. The physics behind this inversion can 
be understood in terms of the different opacities for optical and infrared photons. The 
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opacity of the impinging stellar energy (optical) can be significantly different then the 
opacity of the reprocessed infrared radiation from the planet. Although a one-dimensional 
wavelength dependent radiative transfer model implicitly accounts for these effects by 
modelling the impinging stellar spectrum, full radiative hydrodynamical models must 
use average grey opacities due to computational considerations. In order to encapsulate 
this physics within a full hydrodynamical model I have calculated separate grey opac- 
ities for the incident stellar light and the reprocessed light. Both are still functions of 
the local temperature and pressure, but the stellar opacities are weighted toward optical 
wavelengths. 

Furthermore, to allow for the inversion I have de-coupled the thermal and radiative 
energy components. Each can now evolve independently, though they are linked by the 
appropriate emission and absorption terms. Although computationally more intensive, 
this is a significant improvement to previously employed models. The thermal energy 
equation is given by 

^ + (u-V)e =-PV-u-pK P (T,P)[B{T)-cE B ]+ P K A (T,P)F+e- T % (4.1) 

while the independent radiation energy equation is given by 

— - + V • F = P k p (T, P) [B (T) - cE R ] . (4.2) 

at 

As before, e = c v pT is the internal energy density, T is the local gas temperature, and 
c v is the specific heat. In addition, B (T) = 4<tT 4 , ka(T,P) is the absorption opacity 
averaged with the stellar temperature, and F* is the impinging stellar flux. The term 
proportional to B (T) — cEr, represents the exchange of energy between the thermal and 
radiative components through the emission and absorption of the lower energy photons. 
The pKA (T, P) F±e~ T * term represents the higher energy stellar photons absorbed by the 
gas. Given the relatively low local temperatures, I do not consider the energy re-emitted 
into the higher energy waveband. The radiative flux is treated using the flux-limited 
diffusion method described in Equation fj2.4[) . Equations (|4.1|) and (|4.2|) contain three 
separate opacities: the optical absorption opacity ka, the infrared Planck opacity up, 
and the infrared Rosseland opacity kr that appears through the flux term. 

The effect of decoupling the energy equations and including 3 separate opacities can 
best be understood by considering the radiative behavior of Equations (|4.ip and (14. 2[) . 
ie. in the absence of dynamics. If we express the ratio of the stellar absorption opacity 
to the local opacity as (Hubeny et al. (2003)) 

7^^^, (4.3) 

kp (T, P) 

then the steady-state, radial temperature profile can be expressed as 

TV) = TWe- _ Ml , (4.4) 

4 a 

The temperature in the upper atmosphere, near the stellar photosphere (t* ~ 1) will scale 
as jTq, where To is the equilibrium stellar temperature at the planetary semi-major axis. 
The presence of an inversions thus becomes dependent on the radial structure of 7. In 
Figure ([3]) I show the solution to Equations (|4.1[) and (|4.2p in the absence of dynamics. As 
expected, the inversion is high in the atmosphere, where 7 is around or above unity and t* 
is small. Figure ^ clearly demonstrates how the thermal and radiative energies become 
de-coupled, requiring solving the equations separately. Full dynamical simulations are 
not presented here, but will be in explored in detail in Dobbs-Dixon & Lin (2008). 
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Figure 3. The temperature in HD 209458b concentrating only on the radiative solution to 
Equations (|4.ip and (|4.2p . The left hand panel shows the temperature-pressure profile at the 
subsolar point. The inversion in T (solid line) is quite clear. The dashed line shows the decoupled 
radiative energy Er. The orange line shows the temperature profile at the anti-stellar point. 
The lower left panel shows the behavior of 7 with pressure. The right panel shows the entire 
equatorial slice of the same simulation with the star to the right. 



5. Discussion 

In this proceeding I have reviewed results from full three-dimensional radiative hydro- 
dynamical simulations of the atmospheres of close-in planets. Simulations for a range 
of atmospheric opacities and semi-major axis were performed. The results indicate that 
both stellar proximity and opacity are both crucial factors in determining the behavior of 
planetary atmospheres. Increasing planetary opacities yields less efficient advection and 
colder night-side temperatures. Furthermore, planets situated at 1 and 3 day orbital pe- 
riods exhibit a fairly stable banded jet structure, with jet width increasing with rotation 
period. On the other-hand, dynamics of planets with 7 day orbital periods can be much 
more variable, with wandering jets and the formation of vortices on the night-side. 

The decoupling of the thermal and radiative energies, as well as the introduction of 
multiple opacities relevant for impinging optical stellar photons vrs local re-radiated 
infrared photons allows us to explore in greater detail temperature structure near the 
top of the atmosphere. Radiative solutions presented in Section (4) produce temperature 
inversions similar to those observed and demonstrate the feasibility of simultaneously 
studying these radiative features and the resulting hydrodynamics. The existence of an 
inversion depends critically on the ratio of the optical to infrared opacities and the 
behavior of this ratio with radius. 

Given the importance of both optical and infrared opacities it is fruitful to consider 
how we expect atmospheric properties to vary as we change these two quantities indepen- 
dently. Fortney et al. (2008) suggested a classification planets based on the presence or 
absence of TiVO in the upper atmosphere. Although the presence of TiVO in the upper 
atmosphere is controversial, it would provide the necessary absorption opacity to cause 
an inversion. Fortney's classifications of pM and pL classes correspond to those with and 
without detected inversions. Here extend this classification in terms of the individual 
absorption (ka {T, P)) and Planck [up (T, P))opacities explored in this proceeding. 

As explained in Section (3), the night-side tenrperature is governed by the efficiency of 
advection, which in turn is determined by the cooling timescale where the incident radi- 
ation is deposited. The cooling timescale is directly related to the local Planck opacity 
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Figure 4. The gross behavior of irradiated atmospheres as you change the absorption and local 
Planck opacities. Increasing optical absorption opacity (ka(T, P)) is given in the horizontal 
axis, while the vertical axis gives the local opacity to re-radiated photons (sp (T,P)). The 
diagonal line demarcates the division between planets with an inversion (7 > 1) and those 
without (7 < 1), where 7* = ka (T,P) /up (T,P). Phase measurements coupled to inversion 
determinations will uniquely place you in a given quadrant. 

Kp (T, P) ; large Kp (T, P) will result in a large day- night temperature differential. Such 
a distinction is observable from current and ongoing phase monitoring programs. On the 
other hand, the absorption opacity for the incident stellar light ka {T, P) will play an 
important role in the presence or absence of a temperature inversion. As described in 
Section (4) the ratio 7 = C^ A ( TP l ) 1 ^ 4 will determine the temperature structure in the 
upper atmosphere. If the radial structure of r ye~ T * is such that it peaks in the upper at- 
mosphere there will be an inversion. In these terms, the presence of TiVO (a pL planet) 
would lead to large ka (T, P) , but an inversion would occur only if Kp (T, P) does not 
also increase. In Figure (|4|) I have unified these two concepts onto a single diagram. A 
measurement of both the magnitude of day-night temperature differential and a determi- 
nation of the presence or absence of an inversion will uniquely tell you both the optical 
and infrared opacities and place the planet in one of the quadrants. Such a determination 
will ultimately tell us much about the formation and evolution of these planets. 
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